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Molecular dynamics simulations of 50 Fe grain boundaries were used to understand their interac- 
tion with vacancies and self-interstitial atoms at all atomic positions within 20 A of the boundary, 
which is important for designing radiation-resistant polycrystalline materials. Site-to-site variation 
within the boundary of both vacancy and self-interstitial formation energies is substantial, with the 
majority of sites having lower formation energies than in the bulk. Comparing the vacancy and 
self-interstitial atom binding energies for each site shows that there is an energetic driving force for 
interstitials to preferentially bind to grain boundary sites over vacancies. Furthermore, these results 
provide a valuable dataset for quantifying uncertainty bounds for various grain boundary types at 
the nanoscale, which can be propagated to higher scale simulations of microstructure evolution. 

PACS numbers: 61.72.jd, 61.72.jj, 61.72.Mm, 61.82.Bg, 68.35.Dv 



Defects (vacancies and self-interstitial atoms) caused 
by radiation damage play an important role in the degra- 
dation of physical and mechanical properties of materials 
for nuclear applications [1 j. Future design of nuclear ma- 
terials requires understanding the basic nature of radia- 
tion damage and its evolution in materials at the atomic 
level. Radiation damage begins with the transfer of ki- 
netic energy to a single primary knock-on atom, which 
then collides and displaces neighboring atoms, creating 
an atomic cascade of atoms that are displaced from their 
equilibrium lattice positions. Following this initial event, 
equal numbers of vacancies and self-interstitial atoms 
are generated as atoms attempt to return to an equi- 
librium position within the lattice. Since vacancies, self- 
interstitial atoms, and their clusters have such a profound 
effect on the physical and mechanical properties of these 
alloys, understanding how defects form, diffuse, and are 
absorbed by sinks is crucial to understanding radiation 
damage in fuel cladding materials. Since fuel cladding 
materials are polycrystalline in nature, grain boundaries 
are a significant sink source for defects. 

Molecular dynamics (MD) simulations in single crys- 
tal, bicrystal, and nanocrystal structures have provided 
a fundamental understanding of cascade events and the 
production of defects in polycrystalline materials [2j|3]. 
Atomistic simulations have also been used to examine the 
interaction of defects with grain boundaries [3H9] . For ex- 
ample, MD studies in nanocrystalline metals have shown 
that grain boundaries act as a sink for interstitials [4]. 
While these studies are instrumental to understanding 
the operating mechanisms induced by radiation, these 
types of studies have numerous grain boundaries with 
many confounding effects. For instance, the microstruc- 



ture will contain complex grain boundaries with different 
tilt and twist misorientations, triple junctions, and a dis- 
tribution of grain sizes. It is difficult to ascertain how 
GB character influences the interaction between defects 
and grain boundaries. Therefore, bicrystal simulations 
have been used to explore the influence of grain bound- 
ary (GB) character on specific properties. For instance, 
recently Bai et al. [5 used MD simulations and tempera- 
ture accelerated dynamics to show that grain boundaries 
can act as an effective sink for vacancies and interstitials 
through various mechanisms in Cu. For one such mech- 
anism, interstitials are loaded into the boundary, which 
then acts as a source, emitting interstitials to annihilate 
vacancies in the bulk. While these studies are instrumen- 
tal to understanding the long time scale mechanisms, it 
would be beneficial to use bicrystal MD simulations to 
also understand how GB character may impact the effec- 
tiveness of grain boundaries to act as sinks for defects. 
Additionally, the resulting datasets can be used to pro- 
vide the lower and upper bounds for uncertainty analysis 
in micro- and meso-scale simulations of microstructure 
evolution. 

In this paper, we have investigated the formation en- 
ergies of vacancies and self-interstitial atoms (SIAs) as a 
function of location within/around the GB for 50 (100) 
symmetric tilt grain boundaries in bcc Fe to energetically 
assess the GB sink strength. Nanoscale simulations were 
required to capture the physics of vacancy and intersti- 
tial formation energies at the GB interface. A parallel 
MD code, LAMMPS HQ], was used to run all simula- 
tions in this work. First, a GB database consisting of 
50 (100) symmetric tilt grain boundaries was generated 
using bicrystal simulation cells with 3D periodic bound- 
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ary conditions [cf . [TTHT3] . A minimum distance of 12 nm 
between the two grain boundaries was used to eliminate 
any potential interaction between the two boundaries. 
As with previous work p~2j [13] , multiple initial configu- 
rations with different in-plane rigid body translations and 
an atom deletion criterion were used to properly access 
an optimal minimum energy GB structure via the Polak- 
Ribiere conjugate gradient energy minimization. For an 
initial generation of the structures, the updated version 
of the Mendelev et al. [14] interatomic potential for Fe 
was used. This embedded-atom method [15 potential 
has been shown to perform well in nanoscale simulations 
for nuclear applications [T6] , 

A large number of grain boundaries were used to sam- 
ple the range of GB structures and energies that might be 
observed in polycrystalline materials. The (100) symmet- 
ric tilt grain boundary (STGB) system chosen has several 
low order coincident site lattice (CSL) grain boundaries 
(the £5 and £13 boundaries) as well as both general high 
angle boundaries and low angle grain boundaries (< 15°). 
Figure |TJl shows the GB energy as a function of misorien- 
tation angle for the (100) symmetric tilt grain boundary 
system, similar to that found previously in Fe-Cr sim- 
ulations [17]. The low-order CSL grain boundaries are 
also shown on this figure. For the (100) tilt axis, only 
minor cusps were observed in the energy relationship, 
most noticeably at the £5(310) boundary. In addition 
to many general high angle boundaries, several low angle 
boundaries (< 15°) were also generated. The range of 
GB energies sampled was ~ 500 mJ/m 2 . 

The GB structure plays an important role on the GB 
properties [18]. For low angle boundaries, the grain 
boundary is best represented by an array of discrete dis- 
locations spaced a certain distance apart. However, at 
higher misorientation angles the spacing between dis- 
locations is small enough that dislocation cores over- 
lap and dislocations rearrange to minimize the energy 
of the boundary. The resulting GB structures are often 
characterized by structural units [19]. Grain boundaries 
with certain misorientation angles (and typically a low 
£ value) correspond to 'favored' structural units, while 
all other boundaries are characterized by structural units 
from the two neighboring favored boundaries. Figure 
shows an example for the (100) STGB system, where the 
two £5 boundaries are favored STGBs and the £29 (730) 
boundary is a combination of structural units from the 
two £5 boundaries. The structural units for the £5{210} 
and £5{310} STGBs are labeled B and C, respectively, in 
a naming convention similar to that used for FCC metals 
[TT] . Also notice that the ratio of structural units in the 
£29 GB can be determined by the crystallographic rela- 
tionship of the two favored boundaries. All boundaries 
between the favored structural units of the £5 bound- 
aries and the 'structural units' of the 0° single crystals 
have a similar makeup. 

The formation energies of vacancies and self-interstitial 




15 30 45 60 75 90 

Misorientation angle (degrees) 

(a) 

<100> Symmetric Tilt Grain Boundary Structure 




25(310) 
|B.B| 



229(730) = 2[25(210)] + 1 [15(3 10)] 25(210) 
|BC.BC| l C l 



(b) 



FIG. 1: (a) (100) symmetric tilt grain boundary energy as 
a function of misorientation angle. The low £ grain bound- 
aries in each system are identified, (b) (100) symmetric tilt 
grain boundary structures with structural units outlined for 
the £5(210), £29(730), and £5(310) STGBs. Black and white 
denote atoms on different {100} planes. The different struc- 
tural units are labeled B and C. 



atoms at specific locations within/around the grain 
boundary were calculated for each grain boundary. The 
approach here is similar to that used for modeling point 
defects in Cu [2T] , First, all atoms within 20 A of each 
grain boundary were identified as potential sites. Then, 
for vacancy formation energy simulations, an atom at a 
particular site a was removed and the simulation cell was 
relaxed through an energy minimization. The vacancy 
formation energy for that particular site a was then cal- 
culated by, 



Ef 



^GB 



Eqb + E ( 



coh 



(i) 



where E co h is the cohesive energy /atom of a perfect BCC 
lattice, and Eq B and Eqb are the total energies of the 
simulation cell with and without the vacancy. On the 
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FIG. 2: The vacancy and self-interstitial atom formation en- 
ergies as a function of location for the £5(210), £29(730), and 
£5(310) STGBs. 



other hand, for self-interstitial formation energy simula- 
tions, an atom was inserted approximately 0.5 A away at 
site a, and the simulation cell was again relaxed through 
an energy minimization. The self-interstitial atom for- 
mation energy for site a was calculated by, 



% ^GB 
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(2) 



The only difference between Eq. [T]and Eq. [2] is that E co h 
is subtracted, instead of added, to reflect the added en- 
ergy due to placing that atom in an interstitial site. This 
procedure was then performed for all 50 grain boundaries 
shown in Figure [TJa) and the vacancy and interstitial for- 
mation energies were calculated as a function of location 
for each boundary. 

Figure [2] shows the vacancy and interstitial formation 
energies that correspond to atomic positions in the three 
GB structures shown in Fig. [T}d. The formation energy 
depicted for each location corresponds to a simulation. In 
this graph, the colorbar is normalized such that the high 
value corresponds to the formation energy in the bulk, so 
that vacancy and interstitial energies can be easily com- 
pared. For all three grain boundaries in Fig. [2j there 
are atoms lying along the GB plane that have higher va- 
cancy formation energies similar to the bulk (Ef — 1.72 
eV). However, there are several sites in each grain bound- 
ary that have lower vacancy formation energies than the 
bulk, suggesting an energetic driving force for vacancy 
diffusion from the single crystal to the grain boundary. 
As the vacancy site is shifted away from the boundary, 
the formation energies approach that of the bulk. For in- 
terstitials, most GB sites depicted here have a lower for- 
mation energy than in the bulk (Ef = 3.52 eV). AtomEye 
is used to visualize the simulation results [20] . 

The evolution of the (a) vacancy and (b) interstitial 
formation energies as a function of distance from the 
grain boundary is shown in Fig. [3] for all 50 grain bound- 
aries. Interestingly, most of the formation energies that 
differ from the bulk occur within 5-8 A from the bound- 
ary center. The majority of formation energies within 
this GB-affected region are lower than in the bulk. How- 
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FIG. 3: Evolution of (a) vacancy and (b) interstitial forma- 
tion energies as a function of distance from the grain bound- 
ary for all 50 (100) symmetric tilt grain boundaries (STGBs). 
Low and high angle boundaries are colored differently. The 
inset image is an example of a low angle boundary. 



ever, for both vacancy and self-interstitial atoms, there 
are locations with higher formation energies than in the 
bulk, with the highest vacancy formation energies near 
the boundary center. Additionally, this plot separates 
the high angle boundaries from the low angle bound- 
aries, which is motivated by previous studies [e.g., [8]. 
For interstitial formation energies, most high angle grain 
boundaries have much lower formation energies than that 
in the bulk. The majority of interstitial formation ener- 
gies near bulk values within 5 A of the boundary occur 
in the single crystal regions between dislocations for low 
angle boundaries. For example, the inset image shows 
an example of a low angle boundary. The visible dislo- 
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cations in the low angle boundary have a local effect on 
formation energies and increasing the dislocation spacing 
(i.e., lower misorientation angle) merely results in shift- 
ing these localized regions further apart. For low angle 
boundaries, the formation energies trend to that of an 
isolated dislocation in a single crystal. 

The grain boundary binding energy for vacancies and 
self-interstitial atoms are plotted against each other for 
each site in Fig. [4] The GB binding energy for a par- 
ticular site a is calculated by subtracting the forma- 
tion energy from the bulk formation energy, E^ inding = 
Ef,buik — Ef- The line delineates sites where interstitial 
binding energy is greater than vacancy binding energies 
(above the line). The large amount of binding energies 
above this line indicates that the system energy is de- 
creased more through interstitials occupying GB sites, 
rather than vacancies. Hence, this letter shows that there 
is an energetic driving force for interstitials to segregate 
to GB sites over vacancies. This study also supports the 
interstitial loading-unloading mechanism [5] by showing 
that it is energetically favorable for interstitials to ini- 
tially 'load' grain boundaries for a wide range of grain 
boundary types. This is significant for nuclear applica- 
tions where radiation damage generates lattice defects 
and grain boundaries act as sinks for both vacancies and 
interstitial atoms. 
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FIG. 4: Grain boundary binding energies for vacancies and 
self-interstitial atoms (SIA) for all sampled sites in this study. 
The line delineates sites where interstitial binding energy is 
greater than vacancy binding energies (above the line). 

In summary, this letter shows that self-interstitial 
atoms have a larger energetic driving force for binding 
to grain boundaries than vacancies, based purely on for- 
mation energies of these point defects in thousands of 
GB sites. Fifty (100) symmetric tilt grain boundaries 
were utilized here to sample some of the variability in 



GB degrees of freedom that is observed in experimental 
polycrystalline materials, including general low and high 
angle grain boundaries as well as some low- 1! boundaries. 
By iteratively calculating the formation energies of both 
vacancies and self-interstitial atoms at every site within 
20 A of the boundary, this study was able to methodically 
examine the binding energy associated with GB segrega- 
tion of point defects and calculate the energetic driving 
force for point defects to grain boundaries in Fe. 
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